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In this paper we consider the problem of determining the maximum 
fundamental power in a nonlinear capacitance diode, when the charge 
waveform has a given periodicity and (i) varies between prescribed maxi- 
mum and minimum values, (ii) has a presc?-ibed maximum and a pre- 
scribed maximum slope. Under (i) the maximum obtainable fundamental 
power is first determined. The charge waveform is then further restricted to 
contain no higher than second harmonics, so that the diode is being used as 
a frequency doubter, and, the maximum power transfer is determined. The 
maximum power transfer is also determined under (ii). Particular diodes 
considered, are the abrupt-junction and the graded-junclion ones, with oper- 
ation in the forward conduction region being permitted. 

i. engineer's summary 

This section of the paper is a summary which stresses some of the 
contents of the introduction and summary that follow. It is hoped that 
this will make it easier for the engineer who is involved in parametric 
amplifier and varactor design to deduce the relevant applications of the 
results contained in this paper. 

In the first instance it should be emphasized that an idealized problem, 
based on a mathematical model, is considered. The nonlinear capacitor 
is assumed to be isolated from any external circuits, and we do not dis- 
cuss how the power is fed into or taken from the device. Clearly there 
will be some power lost in the external circuit, and the maximum ob- 
tainable fundamental power determined in this paper is only a theoretical 
maximum, but it would seem to be worthwhile to understand this 
theoretical maximum. When the maximum power transfer from the 
first to the second harmonic is considered, the charge waveform, and 
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hence the current, giving this maximum is determined. Clearly there is 
some relative phase between the first and second harmonics in the cur- 
rent, and the reactance of the output circuit must be adjusted so as to 
obtain this relative phase. 

It is also important to stress that some of the results obtained hold 
for a general, i.e., arbitrary single-valued, voltage-charge relationship, 
and are accordingly applicable to any particular such voltage-charge 
relationship in which the engineer may be interested. We have, for 
simplicity, considered just the abrupt-junction and the graded-j unction 
diodes as special cases, and have idealized the voltage-charge relation- 
ship in the forward conduction region, but other particular diodes can 
be considered as special cases of the general results. We discuss below 
the results which are pertinent to the general voltage-charge relationship. 

Firstly, we have derived the functional form of the charge, waveform 
(of given periodicity and varying between prescribed values) which gives 
the maximum power in the fundamental. The charge waveform is 
composed (see (33) below) of intervals in which it takes on either the 
maximum or minimum prescribed value, or else follows a certain curve. 
The form of the curve depends on the voltage-charge relationship and 
involves parameters which are functionals of the charge waveform 
throughout the entire period, and hence are not known a priori. These 
parameters have to be determined for each particular voltage-charge 
relationship, by solving simultaneous transcendental equations. It is also 
necessary to allow for finite jumps in the charge waveform, and (36) 
below must hold at such a jump. Of course, a jump is not physically 
realizable, since it would correspond to an infinite current, and this 
makes it evident that the maximum is a theoretical one, quite apart 
from losses in the external circuit. It does, however, provide an upper 
bound on the maximum realizable fundamental power. 

In view of the fact that the maximum fundamental power has to be 
determined separately for each specific diode, we derive upper and lower 
bounds for the maximum fundamental power, (11) to (13), which apply 
to a general voltage-charge relationship. For a wide class, the ratio of 
the upper to the lower bound is 1.54. It turns out that, for the particular 
diodes considered, the lower bound is quite close to the actual value. 
Further use is made of the charge waveform giving this lower bound, 
when the power transfer from the fundamental to the second harmonic 
is considered, subject to the charge waveform containing no higher than 
second harmonics. A good approximation to the maximum power 
transfer is obtained by taking the Fourier approximation, up to second 
harmonics, and suitably normalizing so that the approximating charge 
waveform has the prescribed maximum and minimum values. 
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In connection with maximizing the power transfer from the funda- 
mental to the second harmonic, we consider the diode to be a harmonic 
generator, there being input power in the fundamental only. In order to 
make the mathematical problem more tractable, it is supposed that the 
entire output is in the second harmonic. Equations (18) and (19) simply 
state that the maximum power output in the second harmonic, when 
there is input power in the fundamental only, is not greater than the 
maximum obtainable fundamental power without such restrictions, 
and is not less than the maximum fundamental power when there is no 
output or input power in the third and higher harmonics. It is assumed 
here that the charge waveform is continuous. We have already discussed 
the maximum obtainable fundamental power. 

The problem of determining the maximum fundamental power when 
there is no output or input power in the third and higher harmonics is 
still not very tractable, without additional restrictions on the charge 
waveform, and it is thus further supposed that the charge waveform 
contains no higher than second harmonics. The maximum subject to 
this additional restriction is obviously not greater than the maximum 
without it. The significant point about this restriction is that there is 
then no power output or input in the third and higher harmonics, what- 
ever the voltage-charge relationship. We thus determine a canonical 
representation of the charge waveform which contains no higher than 
second harmonics and has prescribed maximum and minimum values. 
By suitable choice of the time origin, this representation contains just 
two parameters which lie in a bounded region. 

Now, it is a straightforward matter to compute numerically the funda- 
mental power for any given voltage-charge relationship and a given 
charge waveform. The numerical maximization of this power with re- 
spect to the two parameters in the above canonical representation is 
also a straightforward process. Thus it is clear that the above procedure 
has general applicability. We add that in the numerical maximization 
process, the two parameters which give the approximating charge wave- 
form (obtained from the charge waveform giving the good lower bound 
to the maximum obtainable fundamental power) are used for starting 
values. 

Consideration is also given to the current-limited diode, in which the 
charge waveform has a prescribed maximum value and a prescribed 
maximum slope (corresponding to maximum current magnitude). 
Again, we determine a two-parameter canonical representation for the 
charge waveform containing no higher than second harmonics, and the 
numerical maximization of the fundamental power, for any given volt- 
age-charge relationship, proceeds along the same lines as in the previous 
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case, except that we no longer have predetermined starting values for 
the two parameters. Lack of space has prevented inclusion of the de- 
termination of the functional form of the charge waveform which gives 
the maximum obtainable fundamental power (without restriction on 
the harmonic content of the charge waveform) in the current-limited 
case. 

II. INTRODUCTION AND SUMMARY 

2.1 Introduction 

We will be concerned with various nonlinear capacitance diodes, these 
being characterized by a nonlinear voltage-charge relationship. Specific 
examples are the abrupt-junction diode and the graded-j unction diode, 
which are composed of diffused p-n junctions. In the former case the 
voltage difference, v, across the diode is proportional to the square of 
the stored charge (per unit area), q, i.e., v a q 2 , while in the latter case 
v cc q\ provided, in both cases, that q ^ 0, which implies that operation 
of the diode does not take place in the forward conduction region. Now 
as electric field strength and barrier width increase, creation of electron- 
hole pairs through secondary impact ionization by both holes and elec- 
trons leads to avalanche multiplication, resulting finally in an effectively 
infinite increase of current with added applied voltage, and this is termed 
reverse breakdown. There is thus a maximum voltage v max , and a cor- 
responding maximum value q mnx of the charge density (which may be 
related to v max . through the actual voltage-charge relationship), above 
which it is not desirable to operate the diode. 

We define the normalized voltage V and the normalized charge Q by 

7 --Li Q--i-. (1) 

"max Qmax 

Hence the normalized voltage-charge relationships for the abrupt-junc- 
tion and graded-j unction diodes, operated in the region between forward 
conduction and reverse breakdown, are 



V = 



Q\ (abrupt) 
Qi, (graded) 



£ Q £ 1. (2) 



It is also possible to operate the diodes partially in the forward conduc- 
tion region, corresponding to Q < 0. The voltage is not very dependent 
on the charge in this region and as an idealization we may assume that 
it is zero throughout. A physical restriction is placed on the maximum 
possible current magnitude, in that the electron velocity is limited by 
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lattice scattering. Throughout most of our analysis we replace this 
condition by a limitation on the minimum charge, so that 

q ^ ffmin = -m(q max ). (3) 

Thus, in the forward conduction region, 

7 = 0, -m ^ Q ^ 0. (4) 

We do, however, give some consideration to the current-limited diode 
in which, instead of (3), 

I * I = *max • (5) 

We will consider charge waveforms that are periodic in time, t, with 
angular frequency w. We define the normalized time x and the normalized 
current I by 

x = ut; I = -J— . (6) 

wgWx 

Thus Q(x) is periodic in x with period 2ir and, since i = dq/dt, 

I = f? = Q'(x). (7) 

ax 

The average real and reactive powers (per unit area) in the nth har- 
monic, p„ and r„ , are given by 

i / \ 2 / r 2rlu \ / r*'i<* 



We define the normalized real and reactive powers in the nth harmonic, 
P„ and R„ , by 

(P. + A) - 2t ' ( P- + ^ . (9) 

wg mttJC y mnx 

We will be concerned with the maximization of the real fundamental 
power, under various conditions, and summarize the results below. We 
note that P„ is not affected by a time shift in the charge waveform, but 
it is reversed in sign by a time reversal of the waveform. 

2.2 The Maximum Obtainable Fundamental Power, When the Charge 
Waveform is Subject to Bounded Variation 

The functional form of the charge waveform which, subject to the 
restriction — m ^ Q(x) ^ 1, maximizes the fundamental power, Pi , is 
found for the general voltage-charge relationship, V = V(Q). The 
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specific form is determined for diodes of interest and the corresponding 
value of max Pi , the maximum obtainable fundamental power, calcu- 
lated. Thus, for the abrupt-junction diode operated in the region between 
forward conduction and reverse breakdown, (2), max P x = 0.1387 and 
the charge waveform Q(x) giving rise to this value is depicted in Fig. 1. 
The corresponding value of the reactive fundamental power is Ri = 
2.43. For the graded-j unction diode, operated in the region between 
forward conduction and reverse breakdown, it is found that max Pi = 
0.408, with Ri = 2.48. The charge waveform giving rise to these values 
is depicted in Fig. 2. The abrupt-junction diode is also considered when 
the region of operation includes forward conduction. Thus, from (2) 
and (4), V(Q) = [max(0,Q)f, -m = Q(x) = 1. Fig. 4 depicts max P, 
and the corresponding Ri as functions of m. The charge waveform Q(x) 
which gives these values when m = 1 is shown in Fig. 5. The somewhat 
idealized voltage-charge relationship given by V(Q) = max(0,Q), 
— »i = Q(x) ^ 1, m > 0, may be treated analytically. It is found in 
this case that 

maxP^^m; ft - 2 (m + 2). (10) 

The charge waveform giving these values is composed of Q(x) = 1,0 
and — m in consecutive intervals of x of length 27r/3. 

It is observed that the charge waveform which gives rise to max Pi , 
for the various diodes, contains at least one discontinuity (or jump) in 
a period. A jump, of course, is not physically realizable, since it would 
correspond to an infinite current, so max P t cannot actually be attained. 

Finally, upper and lower bounds are obtained on the maximum ob- 
tainable fundamental power, max Pi , for the general voltage-charge 
relationship V = V(Q), with —m = Q{x) = 1. Thus, it is shown that 

'-^L = maxP, = 4(1 + m)U, (11) 

where 

L= max [(p-r)F(ff) + (T-ir)7(p) + (ff-p)7(r); : (12) 

— m^ (ir,p,T)^ 1 



and 



U = min { max [Xer - V(<r)\ - min [Xcr - 7(cr)]}. (13) 

X -m<,o^l —m<o<l 



Moreover, it is shown that 

L = (1 + m)U = 2L. (14; 
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The bounds in (14) cannot be improved without restriction on V(Q), 
but if [(p + m)V(l) - (m + l)F(p) + (1 - p)V(-m)] does not 
change sign in — m ^ p ^ 1, then L = (1 + m)U and the ratio of the 
upper to the lower bound in (11) becomes 1.54. The class of voltage- 
charge relationships 

V(Q) = [max (0,Q)]\ -m ^ Q ^ 1; m ^ 0, v £ 1, (15) 

which includes the particular diodes considered, satisfies the above 
condition, and in (his case 

L = ro + (l - i) [(1 + m)vV m9 - x \ (16) 

For the particular cases considered, the lower bound in (11) is fairly 
close to max Pi . 

A lower bound is also obtained, for a general voltage-charge relation- 
ship V = V(Q), with -m ^ Q(x) ^ 1, for P x such that P l + P 2 = 0. 
It is shown that 

max [Pi | P, + P 2 = 0] ^ (1.87)L. (17) 

2.3 The Maximization of the Power Transfer in a Frequency Doubter, 
With Bounded Charge Waveform 

Here we are interested in maximizing the power transfer from the 
fundamental to the second harmonic, when the diode is being used as 
a harmonic generator. Thus there must be input power at the funda- 
mental frequency only, i.e., I\ > and P„ ^ 0, n ^ 2. In order to make 
the problem more tractable we suppose that the entire power output is 
put is in the second harmonic, so that P„ = 0, n ^ 3. It follows that 
Pi + P 2 = 0, provided that the charge waveform is continuous, since 
then 2^"~i P» = 0. We observe that 

max [-P 2 | P„ ^ 0, n ^ 3] 

(18) 
^ max [P, | P„ ^ 0, n ^ 3] ^ max P, , 

and 

max [-P a | P„ gO, n £ 3] ^ max [-P, | P„ = 0, n ^ 3] 

(19) 
= max [P, | P„ = 0, n ^ 3]. 

Even the problem of determining max [P x | P„ = 0, >i ^ 3], that is, 
max Pi subject to P n = 0, n ^ 3, is not very tractable, without addi- 
tional restrictions on the charge waveform. Thus, it is supposed that 
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the charge, and hence the current, contains no higher than second har- 
monics. The conditions Pi + P 2 = and P„ = 0, n ^ 3, are then iden- 
tically satisfied, independently of the voltage-charge relationship V = 

V(Q). 

Now, a change in the time origin does not affect the power transfer. 
Hence, the canonical representation of a charge waveform which con- 
tains no higher than second harmonics and is such that Q(x) = Qmin = 
—m and Q(2 tan -1 s) = Q,„ as = 1, is constructed. In addition to the 
parameter s there is the parameter y which is subject to the restriction 
^ y ^ (1 — s 2 ), which of course also implies that s 2 ^ 1. It is found 
that Pi = on y = and on y = (1 — a 2 ), independently of the vol- 
tage-charge relationship. Moreover, Pi(s,y) = -Pi(—s,y) and in 
particular P y = on s = also, so that it is sufficient to consider only 
the region —1 ^ s ^ 0, ^ y ^ (1 — s 2 ) and to maximize | Pi |. The 
abrupt-junction diode, operated in the region between forward conduc- 
tion and reverse breakdown, may be treated analytically, and it is found 
that the maximum power transfer is 0.281, as compared with the max- 
imum obtainable fundamental power of 0.687. The corresponding re- 
active fundamental powers are 1.46 and 2.43, and the charge waveform 
giving the maximum power transfer is depicted in Fig. 6, which should 
be compared with Fig. 1. 

In order to determine the maximum power transfer for a general 
voltage-charge relationship, recourse must be made to numerical com- 
putation. However, a prior step is the determination of a charge wave- 
form which provides a reasonable approximation to the maximum power 
transfer, and hence provides starting values for s and y in the numerical 
maximization process. A good lower bound was obtained for the max- 
imum obtainable fundamental power. Furthermore, for a wide class of 
voltage-charge relationships V = V(Q), the charge waveform Q(x) 
giving this lower bound satisfies Q mnx = 1 and Q mi „ = -m. The class 
of voltage-charge relationships (15) falls within this class. Thus it would 
seem feasible that a reasonable approximation to the maximum power 
transfer will be obtained by taking the Fourier approximation, up to 
the second harmonics, of the charge waveform giving the good lower 
bound for the maximum obtainable fundamental power, and suitably 
shifting and expanding (or contracting) the Fourier approximation so 
that the resulting charge waveform Q(x) satisfies Q max = 1 and Q, llin = 
— m. This is the procedure adopted and, for the abrupt-junction diode, 
operated in the region between forward conduction and reverse break- 
down, it actually yields the charge waveform that gives the maximum 
power transfer. 
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The results of the numerical maximization process are tabulated in 
Section 5.4. Tables II and III, for the cases v = 2 and v = f , in (15), show 
the values of max Pj , the maximum power transfer, and the correspond- 
ing values of R, and R 2 , the reactive powers in the fundamental and 
second harmonic, I max , the maximum normalized current magnitude, 
and (6 + c") and (d 1 + e )> the squares of the amplitudes of the first 
and second harmonics in the charge waveform, for several values of m. 
Tables IV and V show the values of — s and y which give max l\ and 
also y and iV", the value of P, corresponding to the starting values 
y and — .s ( " = l/\/3. It is interesting to observe how close Pi is 
to max l J i , particularly for the smaller values of m. Table VI compares 
max Pi with the maximum obtainable fundamental power, max P t , in 
the case v = 2, for several values of m. It is also worth noting that in 
the case v = |> m «* we have max Pi = 0.162, whereas max Pi = 
0.408. 

2.4 The Maximization of the Power Transfer in a Frequency Doubler, 
for the Current- Limited Diode 

We finally turn our attention to the current-limited diode in which 
(5), instead of (3), holds. Thus, from (5) to (7), 

| Q'(x) | ^ -^4 = K - . (20) 

For the P N abrupt-junction diode of germanium 1,2 

y max ~ 1.03 X 10 ,3 (iVT ' 726 volts, 

(2D 
i max ~ 1.6 X 10 N amps/cm", 

where JV is the donor concentration in cm -3 . But, from the voltage- 
charge relationship, 

q m J = 2c e Nv max , (22) 

where e denotes electron charge. Hence, 

7milx ~ 2.16 X 10 9 (JV) - ,37S coulombs/cm 2 , (23) 

and 

K = l J^L ~ 0.74 X 10- 3 (AT) - 8625 sec -1 . (24) 

Qm ux 

For JV = 2 X 10 lfi , a reasonable value, k « 10 11 sec -1 , which is in the 

range of angular frequencies of interest. 
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We consider the problem of maximizing the power transfer from the 
fundamental to the second harmonic, when the diode is being used as 
a frequency doubler, and, as previously, the additional assumption is 
made that the charge waveform Q(x), and hence the current, contains 
no higher than second harmonics. The first step is the construction of 
the canonical representation of Q(x) such that Q max = 1, Q'(v) = 
Qmin - — * and Q'(2 tan -1 s) = QLx ^ ft- In addition to the param- 
eter s there is the parameter y which is subject to the restriction ^ 
y ^ i(l — s 2 ), which of course also implies s 2 ^ 1. It is found that 
Pi = on y = |(1 _ s *)> independently of the voltage-charge rela- 
tionship. Since, if Q(x) = Q(ir — x), then Q max = 1, Q max = k and 
QLin ^ — k, it is sufficient to consider the above canonical representa- 
tion and to maximize | Pi |, in order to maximize P, subject to Q mu = 
1, | Q' | ma x = k- We denote this maximum by n(A;). For the abrupt- 
junction diode operated in the region between forward conduction and 
reverse breakdown, the determination of II (k) is carried out analytically 
for k sufficiently small that Q inin ^ 0. It is found that Tl(k) = 0.731& 3 , 
for ^ k ^ 0.(581. Combining this result with that obtained when the 
charge waveform is subject just to bounded variation, ^ Q(x) ^ 1, 
it is shown that, from the viewpoint of maximizing the actual funda- 
mental real power p, , the optimum operating frequency lies in the range 

j.299 ^ ( " 9n "' x) ^ 1.468, (25) 

'max 

and that 



54 (max p t ) < 2(4)' 1 QQ 



1 ^ ^ mttA w ^ *_y^j_ < 106 (2fi) 

l, 'max'' max ^ O 

For the abrupt-junction diode which is allowed to operate partly in the 
forward conduction region, the maximization of the power transfer is de- 
termined by numerical computation. For the values of s and y which give 
max | Pi |, i.e., Tl(k), the reactive powers Ri and R 2 , and Q min , i.e., 
— M(k), were calculated, the results being given in Table VII (Section 
6.4). It is shown that max Pi subject to Q max S 1 and | Q' \ mox ^ k is at- 
tained with Qn.nx = 1 and | Q' | mMt = k. For k < 0.681 it can also be at- 
tained with 1.468ft ^ Q,„ax < 1 and | Q' | m ax = /v- Optimizing with respect 
to the frequency it appears that 20 (max pi) ~ i ma xV ma x . Thus a consider- 
able improvement is obtained by permitting operation in the forward 
conduction region. The optimum frequency in this case is roughly one- 
fifth that in the case when operation is not allowed in the forward con- 
duction region, although close to max pi may be obtained at one-third 
the frequency. 
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In conclusion, we add that lack of space has necessitated the omission 
of several aspects of this problem, and in particular of the determination 
of the maximum obtainable fundamental power when the periodic 
charge waveform is restricted only to have bounded slope. 

III. THE CHARGE WAVEFORM WHICH, SUBJECT TO BOUNDED VARIATION, 
MAXIMIZES THE POWER IN THE FUNDAMENTAL HARMONIC 

3.1 The Functional Form of the Charge Waveform 
From(l), (0), (7), (8) and (9), 

P« + jR n = (T Q\x) «-*" dx\ (T V[Q(x)] e ini dx\ . (27) 

It is noted that P n is not affected by a time shift in the charge waveform 
Q(x), but it is reversed in sign by a time reversal of the waveform. On 
the other hand, R„ is not affected by either a time shift or a time re- 
versal in the charge waveform. Integrating by parts the first integral 
in (27), and remembering that Q(x) is periodic with period 2w, and 
then separating real and imaginary parts, 

P n = n(a„S„ - 0„7n); Rn = n(a„y n + 0„8 n ), (28) 

where 



a„ = 

'0 



/ Q(x) sin nx dx; 0„ = / Q(x) cos nx dx; 

r'2x />2ir 

In = / V^[Q(aO] sin nx dx; 8„ = / V[Q(x)] cos nx dx. 
From (28) and (20) we may express P„ as a double integral, 

1 P»= I [ Q(x)V[Q(y)] sin n(x - y) dx dy. 

n Jo Jo 



(29) 



[30) 



To find the functional form of Q(x) which, subject to the restriction 

-m ^ Q(x) SI, (31) 

maximizes P { , we set 

Q(x) = \{\ + m) -sech R(x) - m\, (32) 

so that the inequalities in (31 ) are satisfied. A variational procedure 
applied to (30) then shows that for stationary values of Pi , we have, 
for each x, 
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Q( x ) = -m 7 or Q(x) = 1, or 

, , ( 71 cos x - 5i sin x) (33) 

7 [Q{X)] = (a^osx-frsmx) ' 

where ay , ft , 71 , and 5! are as defined in (29). This, then, is the func- 
tional form of Q(x) which maximizes Pi . Evaluation of the integrals 
in (29) will lead to four equations for the four unknowns «i , ft , 7i > 
and 5i . Note, however, that 

d R71 cos x — l\ sin x)~\ («i3i — ffm) 

rfx L(«i cosa; - /3i sin.r)J («i cos x - /3, sin .r) 2 



-Pi 



(34) 



(ai cos x — ft sin .t) 2 ' 



from (28), is of one sign. Since we are not interested in T\ = 0, which 
case arises in particular if Q(x) = const, it follows that allowance must 
be made for discontinuities in Q(x), since we require that Q(x) be peri- 
odic. Supposing that Q(x) is discontinuous at x = <p, we obtain a con- 
dition by integrating the equation 

V'[Q(,)]Q'M = ! TI M * ~ * si " *\ Q'M, (35) 

(ai cos x — /3i sin x) 
from x = if - to a; = v» + 0. This gives 

\V[Q(x)\Y,tl = <*"■'«'- ^'j [«■)]« ■ (36) 

L * J * (ai COS *» — |8i Hill ^) 

3.2 77m? Cfeor^e Waveform for the Abrupt-Junction Diode 

In normalized form the voltage-charge relationship for the abrupt- 
junction diode operated in the region between forward conduction and 
reverse breakdown is 

V(Q) = Q\ ^ Q(x) g 1, (37) 

so that m = in (31). We make use of the fact that Pi is invariant 
under the transformation Q(x) •* Q(x - 6), and choose 6 so that ft = 
0, since this leads to a simplification of the analysis. Let us define a and 
b by the equations 

7l = 2aa, , 5! = 2fea, ; ft - 0. (38) 

Then, from (28), 

Pi = 2ba{. (39) 

It is clear that max P t > 0, and hence that b > 0. The functional form 
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of Q(x) for max P, is, from (33), (37) and (38), 

Q(x) = 0, or Q(.v) = 1, or Q(x) = (a - b tan x). (40) 

Rejecting combinations which lead to Pi = 0, we are led to the con- 
clusion that, within a cycle, Q(.v) = 1 for an interval, it then follows 
the curve Q(x) = (a — b tan x) and then Q(x) = for an interval, 
after which it jumps from to 1 and the cycle is repeated. 

Let <p be a value of x at which a jump in Q(x) from to 1 occurs. 
Then (30), (37), and (38) give 

(2a - 1) 



tan ip = 



2b 



[41) 



Thus we obtain max I\ by taking 

1, for v < x ^ tv + tan" 1 [(a - l)/6]; 

a - b tan x), for t + tan -1 [(a - l)/6] ^ x 
^ 7t + tan -1 (a/6) ; 
0, for 7T + tan -1 (a/b) gx < 2tt + p, 



QU) = 



(42) 



where 



and 



-£ < tan"' [(a - l)/ft] < ^ < tan -1 (a/6) < J , 



Q(a: + 27r) = Q(x), all*. 



(43) 



(44) 



Now ati , j8i , 7i , and Si may be calculated from (29), (37), (42) and 
(44). Substitution into (38) then leads to 



(2a - 1) cos*» = 26{[(a - 1 r + 6 2 ] 1 - («'"' + &T) J 
26 cos*> + sin *> + 36t - {(a + l)[(a - l) 2 + 6 2 ] } 

- a(a 2 + 6 2 ) 5 |; 
sin*) = {[(a-l) 2 + 6 2 ]'-(a 2 + 6 2 ) 5 }, 
where 
r = 6[tanh _1 [a(a 2 + 6 2 ) _i ) - tanh -1 { (a - 1 )[(a - l) 2 + 6 2 ]" 3 )]. 



(45) 



(40) 



It would appear that we now have one too many conditions on a, 6 and 
ip because of the relationship in (41 ), which was obtained from the jump 
condition at x = ip, but it is observed that the first and last equations 



cm 
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(47) 



in (45) are consistent with (41). Since \<p\ < tt/2 and b > 0, (41) gives 

cos*. = 2b[(2a - l) 2 + 46 2 ] - *; 

sin^= (2a - l)[(2o- l) 2 + 4b 2 ]"*: 

Substituting into the first equation in (45), we obtain 

(2a - \)[(2a - 1 ) 2 + 4b 2 ]-* = {[(a - l) 2 + b 2 ] 4 - (a 2 + b 2 )*}. (48) 

A solution to (48) is a . = \ and, moreover, this is the only solution 
since if a >\ the L.H.S. > and the R.H.S. < 0, and vice versa. Thus, 

a = i <p = 0. (49) 

The second equation in (45), using the definition of t given in (46), 
now leads to an equation for b, namely 

3b 2 tantT 1 [(1 + 4b 2 )" 1 ] = ft(i + 4b 2 ) 1 - b], (50) 

and (HI)) and the expression for ai give 

Pi = 26{1 + 2b tanh- 2 [(l + 4b 2 )-*]} 2 = A [2b + (1 + 4fc 2 )'l 2 , (51) 

using (50). Equation (50) was solved numerically and it was found that 

b = 0.14136; max Pj = 0.6868. (52) 

The shape of Q(x) which gives this maximum value of l\ is shown in 
Fig. 1. From (28) and (38) the corresponding reactive fundamental 




Fig. 1 — Charge waveform for maximum obtainable fundamental power in 
abrupt-junction diode operated in the region between forward conduction and 
reverse breakdown. 



POWER IN NONLINEAR CAPACITANCE DIODES G91 

power is given by 

R x = 2aa{ = iJL = 2.429. (53) 

Note that the reactive power is about three and a half times as large 
as the real power. 

3.3 The Charge Waveform fur the Graded- J unction Diode 

We now turn our attention to the spcond diode of interest, namely 
the graded-j unction diode, and suppose that it is operated in the region 
between forward conduction and reverse breakdown. In normalized 
form the voltage-charge relationship is 

V(Q) = Q\ S Q(x) ^ 1. (54) 

The determination of the maximum obtainable fundamental power, 
max 7 J ! , is carried out along the same lines as for the abrupt-junction 
diode, although the details are more involved. The analytical form of 
the charge waveform Q(x) which gives max P x is 



Q(x) = 



where 



1, for $ < x ^ ir + tan -1 [(a - \)/b\\ 
(a - b tan x)\ for n + tan" 1 [(a - l)/b] ^ x 

^ tv + tan -1 (a/6); 
0, for 7T + tan -1 (a/6) ^ x < 2w + ^, 



(55) 



-| < tan -1 [(a - l)/6] < + < tan _1 (a/6) < £, (56) 

and (44) holds. Here 

3 3 

7, = = aai ; 5, = - baa ; ft = 0, (57) 

which leads to three equations for a, b and \p. These equations are con- 
sistent with the jump condition (36) which gives 

tan ^ = (3a - 2)/(36). (58) 

Elimination of \J/ leads to two equations for a and 6 which were solved 
numerically, giving 

6 = 0.1101)8; a = 0.G7375. (59) 
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These values lead to 

max P x = 0.4084; R, = 2.479. (60) 

The corresponding charge waveform Q(x) is depicted in Fig. 2. 

3.4 The Abrupt- J unction Diode When the Region of Operation Includes 
Forward Conduction 

In this case the normalized voltage-charge relationship is, from (2) 
and (4), 

V(Q) = [max (0,Q)f, -m fe Q(*) SI, m > 0. (61) 

As previously, we translate Q(.r) so that ft = and again define, a and 
6 by (38), so that (39) for P t also holds. From (33), (38), and (61), 
the functional form of Q(x) for max Pi is 



Q(x) = -m, or Q(x) = 1, 
or max [0,Q(x)\ = (a — It tan x). 



(62) 



Thus we are led to the conclusion that within a cycle Q{x) = 1 for an 
interval, it then follows the curve Q(x) = (a - b tan x) until the point 
at which Q(x) = where it jumps to the value —m, and after Q{x) = 
—m for an interval it jumps to the value 1 and the cycle is repeated. 
Thus in this idealized case there are two discontinuities in Q(x) in one 
cycle. Note that according to (36), together with (38), the jump of 
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Fig. 2 — Charge waveform for maximum obtainable fundamental power in 
graded-junction diode operated in the region between forward conduction and 
reverse breakdown. 
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Q(x) from to — m occurs at x = \p where (a — b tan \f/) = 0, since 
V(0) = and V( — m) = 0. Hence we obtain max P t by taking 

1, for ip < x ^ 7T + tan -1 [(a - l)/6]j 

(a - 6 tan .r), for tt + tan -1 [(a - l)/6] ^ .r 

< 7T + tan -1 (a/6); 

— m, for tt + tan -1 (a/6) < x < 27r + <p, 



Q(x) = 



(03) 



where (43) and (44) hold. In this case the jump condition at. x = <p, 
gives 

[2a (1 + m) - 1] , . 

tan <p = — r . (04) 

26(1 + m) 

The calculation of ai , ft , ?i , and 5i , and substitution into (38), 
leads to three equations for a, 6, and </?, which are consistent with (64). 
The elimination of p leads to two equations for a and 6, which quan- 
tities of course are functions of m. It was found to be possible to elim- 
inate m analytically from these two equations, so that instead of solv- 
ing the two simultaneous equations for a and 6 for given values of m, 
the single relation between a and 6 which did not involve m was solved 
for 6 for given values of a. Thus a parametric solution was obtained in 
the form 6 = 6(a), m = m(a). From this a and 6 were plotted graph- 
ically against m and the results are shown in Fig. 3. It was shown ana- 
lytically that 

1 < 4a(l +m) ^ 2, (05) 

the upper bound being attained for m = and the lower bound being 
approached for m —* » . Also, as m — *■ =o it is found that 

6 ~ V3a; max P, ~ '^^ m; R } ~ ? m, (66) 

where R\ is the reactive power in the fundamental. Fig. (4) shows 
max / J ! and the corresponding R x as functions of m. It is interesting to 
note that the ratio (max Pi)/R\ increases with increasing m from its 
initial value of 0.28, its asymptotic value being -\/3, from (00). The 
charge waveform Q(x) giving rise to max P x is shown, for m = 1, in 
Fig. 5. 

3.5 The Charge Waveform for an Idealized Voltage-Charge Relationship 

We now consider a special voltage-charge relationship which may be 
handled analytically. Thus we suppose that the capacitance has a finite 
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ig. 3 — Parameters in charge waveform for maximum obtainable fundamental 
or in abrupt-junction diode operated partly in forward conduction region, vs. 
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Fig. 5 — Charge waveform for maximum obtainable fundamental power in 
abrupt-junction diode operated partly in forward conduction region. 



constant value for reverse bias and is infinite for forward bias, and 
hence in normalized form 

V{Q) = max (0,Q), -m ^ Q(x) g 1; m > 0. (67) 

Since V'(Q) is constant except possibly at Q = 0, where it is indeter- 
minate, we deduce from (33) that Q(x) has one of the values 1, 0, and 
— m at each point. Omitting further details, it is found that max Pi is 
given by 



1, < x < 2t/3; 

'/' : (0, 2x/3 < x < 4tt/3; 

— TO, 47r/3 < X < 2ir. 



(68) 



Also, 



o 3V3 

max Pi. = — -— m; 



R 



1 ,= |(to + 2). 



(69) 



Note that, as might be expected, these values are asymptotically, as 
m -* « , the same as for the voltage-charge relationship in (61), as is 
seen from (66). 



696 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1962 

IV. BOUNDS ON THE MAXIMUM OBTAINABLE FUNDAMENTAL POWER 

4.1 Lower Bounds 

We now derive some lower bounds for the maximum obtainable power 
in the fundamental, for a general voltage-charge relationship, by the 
simple expedient of choosing specific charge waveforms. Any Pi which 
we obtain is, of course, a lower bound for max Pi . Thus we consider 
the charge waveforms 

IV, on Tr ; 

Q( x ) = < p , on T, ; (70) 

t, on r 3 , 

where each Tj(j = 1,2,3) is a finite collection of nonintersecting inter- 
vals, open at the left and closed at the right, and furthermore 

r,n r* = o, j*k; u r,= (0,2*]. (7i) 

J=l 
From (28), (29), (70) and (71), 

P„ = nL(a,p,r) 



cos nx dx 1 I / sin nx dx 



( / sin nx dx J ( / cos nx dx j , 



(72) 



where 

L(<t,P,t) = [(p - r)V(a) + (r - a)V(p) + ((T - p)V(r)}. (73) 

The significant point here is that we can choose the intervals r x and 
T 2 to make Pi as large as possible, for the waveform class of (70), in- 
dependently of the functional form of the voltage-charge relationship, 
V = V(Q). This is still true if we wish to make P, as large as possible 
subject to the condition Pi + P 2 = 0, say, since the factor containing 
V, namely L(<t,p,t), occurs in each P„ . Note, from (73), that L(<t,p,t) 
vanishes unless a, p, and r are unequal. Also, if (<r,p,r) undergo a cyclic 
permutation then L(<t,p,t) is unaltered, but if (<t,p,t) undergo an anti- 
cyclic permutation then L{a,p,r) is reversed in sign. We suppose that 
the charge waveform has bounded variation as in (31) and define 

L - max [L(<t,p,t)} ^ 0. (74) 



(77: 
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Also from (72) it is seen that Pi changes sign if Ti and r 2 are inter- 
changed, which is equivalent to an anticyclic permutation of (a,p,r). 

Thus wo are interested in making the modulus (or magnitude) of the 
bracketed expression following the factor L(<j,p,t) in (72) as large as 
possible, in order to obtain as large as possible a lower bound for 
max Pi . We will restrict ourselves to special Tj and find the maximum 
modulus of the bracketed expression in (72) for these subclasses. In par- 
ticular, we consider 

r, = (o,\] ; r 2 = (ixM, o < \ ^ n < p < 2t. (75) 

Then, from (72), 

P„ = - U*, P ,T)F{n\,nv,nv) , (76) 

n 

where 

F(\,n,v) = [sin (v — X) — sin (p. — X) + sin p. — sin v\ 

= 4 sin [(v - fi)/2] sin (A/2) sin [(p + p - X)/2]. 

We first set X = /x and determine p and v to maximize F(p,p,v) which 
from (75) and (77) is seen to be positive. The stationary values of 
[sin {v — fi) + sin n — sin v\ arc given by 

cos n = cos (ft — v) = cos v. (78) 

Hence F(n,n,v) is a maximum for p. — 2tt/3, v = 4x/3 and from (74), 
(76), and (77) the corresponding maximum of P x is 

P 1 = ?Vj*L. (79) 

Now for the voltage-charge relationship (37) it is readily verified that 
L(<t,p,t), as defined in (73), has a maximum value of \ which is attained 
for a = 1, p = \, t — 0, and hence in this case we obtain the value 
P\ = 0.650, which is quite close to the value of max Pi given in (52). 
We now consider the maximization of F(\,p,v) subject to the con- 
dition 

F(\ji,v) + iP(2X,2 M » = 0, (80) 

corresponding to P x + P 2 = 0. Using the second part of (77), (80) 
becomes 

1 + 4 cos [( v - n )/2] cos (A/2) cos [(» + /* - A )/2| = 0, (81 ) 
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supposing that F(\,n,v) ^ 0. It is interesting to note that (81) cannot 
be satisfied with X = /*. It is found that F(\,n,p) is maximized, subject 
to (75) and (81), by 

X = 2(tt-0); n = 6; v = (2ir - 0), (82) 

where 

cos d = -(£)», (2t/3 < 6 < tt), (83) 

and the corresponding value of Pi , with Pi + P2 = 0, is 

P,/(4L) - [1 - (*)*]» = 0.468. (84) 

4.2 An Upper Bound, and its Relationship to a Lower Bound 

In Appendix A we give the derivation of an upper bound, for a gen- 
eral voltage-charge relationship, on the maximum obtainable funda- 
mental power, using the fact that the charge waveform is of bounded 
variation, (31). It is shown that 

maxF, ^ 4(1 + m)U, (85) 



where 



U = min( max [X<r - V(a)\ - min [Xt - V(<r)]}. (80) 



In the previous section we showed, by example, that 

maxP.^'^L, (87) 

where L is defined by (73) and (74). From Appendix A, we have 

1 ^ (1 4- m)U/L ^ 2, (88) 

and these bounds cannot be improved without restriction on the voltage- 
charge relationship. However, there is a large class of voltage-charge 
relationships for which the lower bound is attained, namely those for 
which [( P + m)V(l) - (m 4- l)V(p) + (1 - p)V(-m)] does not 
change sign in — m ^ p ^ 1. From (85) and (87) it follows that 

3yg ^maxP! ^ 4> . f L =(i +m )U. (89) 

2 Li 

Also, for the above class, L in (74) is given with a = 1, t = — m, or 
vice versa, and U in (86) is given with X = [V(l) — V(-m)]/(l +m). 
A class of voltage-charge relationships of interest is 

V(Q) = [max (0,Q)}\ -m £ Q £ 1; m ^ 0, v ^ 1, (90) 
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of which we have already considered the cases v = 2, v = 1, and v = 

\ (with m — 0). It is readily seen that this class satisfies the above 
condition, and hence 



L = max {(p + m) — (1 + m) [max (0,p)]"} ; 

-mgpgl 



C7 = max 



-mgpgl 1(1 + m) 



- [max (0,p)]' 



(91) 



— mm 



-mgpgi {(I + m) 



- [max (0,p)]' 



Thus, 



L = w + 



H) 



i)[(l + m)H 



-i/(i—i) 



= (1 +m)C/, (92) 



and the bounds on (max P\)/L in (89) hold. From (69) and (92) it is 
seen that the lower bound is exact for the case v = 1, m > 0. For v = 
f and m = 0, 3\/3 L/2 = 0.385 as compared with max P x = 0.408. 
For v = 2, L = {2m + 1) 2 /[4(1 + w)], and Table I shows the ratio 
2 (max Pi)/(3v / 3L) for several values of m, and it is noted that the 
lower bound improves with increasing m. 







Table 


l-(p = 


2) 






m 

2(max P,) 

3V3L 



1.058 


0.589 
1.042 


1.20 
1.027 


1.89 
1.018 


3.07 
1.012 


5.50 
1.006 



V. THE MAXIMIZATION OF THE POWER TRANSFER FROM THE FUNDA- 
MENTAL TO SECOND HARMONIC, WITH BOUNDED CHARGE WAVEFORM 

5.1 The Canonical Representation of the Charge Waveform 

We wish to consider the problem of maximizing the power transfer 
from the fundamental to the second harmonic, when the charge wave- 
form contains no higher than second harmonics, so that 

Q(x) = a + b sin x + c cos x + d sin 2x + e cos 2.r. (93) 
We also impose the conditions 

Qmnx = 1; Qmin = ~m. (94) 

Note that it does not follow a priori that the maximum power transfer 
subject to (94) is equal to the maximum subject to Q max ^ 1, Qmin ^ 
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— to. We observe, however, that for the voltage-charge relationship 
(90), 

max [Pi | Q,„ ax = r, Q miB = — q\ 

= r (r+1) max [P, | Q max = 1, Q min = -ff/r], < ;• ^ 1, 

as may be deduced from (28) and (29). Thus it is sufficient to deter- 
mine max [Pi | Q max = 1, Qmin = — to], that is, max Pi subject to the 
conditions of (94), for a range of values of m. 

A canonical representation of Q(x) is found in Appendix B. In addi- 
tion to the two conditions in (94) it is supposed, by a suitable choice 
of time origin, that 

Q(t) = Q mm = -m. (96) 

Thus the five coefficients in (93) arc given in terms of two parameters 
and it is found that 

a = [(c — e) — to]; 6 ■» 2d = (1 + m)sij; 

c = (1 +»)[J(1 -s 4 ) -s 2 y]; (97) 

e - |(1 +»)fo(l -s 2 ) - 1(1 4-s 2 ) 2 ]. 

The parameter s arises from the equation 

Q(2tan _1 s) = Q max = 1. (98) 

The parameter y is subject to the condition 

^y ^ (1 - s 2 ), (99) 

which of course also implies that s ^ 1. Thus we have a two-parameter 
canonical representation of Q(x), and these two parameters lie in a 
bounded region. Moreover, it is shown in Appendix B that, independ- 
ently of the voltage-charge relationship V = V(Q), 

P, \ tt=0 = 0; P x L /=(1 _ s2) = 0, (100) 

so that Pi vanishes on the boundary of this region. Also it is seen, from 
(93) and (97), that changing the sign of s is equivalent to the trans- 
formation Q(x) — * Q(2t — x), and hence 

Pi(s,y) = -Pi(s,y); Pi |„=o = 0. (101) 

5.2 The Abrupt- J unction Diode 

We now consider the abrupt-junction diode operated in the region 
between forward conduction and reverse breakdown. From (28), (29), 
(37), (93) and (97), with m = 0, it follows that 

Pi= -js(l+s 2 )V[(l -s 2 ) -y\. (102) 
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The maximum of (102) subject to (99) is 



max Pi = 



4tt 



81\/3 



= = 0.2814, 



(103) 



being given by s = — (l/v3), U = \- Thus a charge waveform giving 
max 1\ is 

Q(x) = \ + ^ [2 sin (x + ^) + sin 2 (x + ^j , (104) 

and (lie corresponding fundamental reactive power is found to be 

R x = ^ = 1.462. (105) 

Q(x) = Q(x — (2t/3)) is depicted in Fig. 6. It is interesting to com- 
pare (52) and (103), and Figs. 1 and 6. We comment that the above 
results may be obtained quite elegantly, without using the canonical 
representation of the charge waveform. 

5.3 A Charge Waveform Which Provides an Approximation to the Maxi- 
mum Power Transfer 

In Section IV wc obtained a lower bound to the maximum obtainable 
fundamental power, (87), and it was seen to be a close bound in the 
particular cases considered. The charge waveform giving this lower 
bound is one which has values a, p, and r on consecutive intervals of 
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Fig. 6 — Charge waveform for maximum power transfer from fundamental to 
second harmonic in abrupt -junction diode operated in the region between forward 
conduction and reverse breakdown. 
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x of length 2t/3. Here a, p, and r are those values which, subject to 
— m ^ (<t,p,t) ^ 1, maximize L{o,p,t), as denned in (73). It was also 
pointed out that if [( P + m)V{l) - (1 + m)V(p) + (1 - p)V(-m)), 
i.e., L(l,p, — m), does not change sign in — m ^ p ^ 1, then L{o,p,r) 
is maximized with a = 1, t = — m (or vice versa) and a suitable value 
of p. The class of voltage-charge relationships given in (90) satisfies 
this condition and then 

p= [(1 + m )„]-*'<-». (100) 

Now the Fourier coefficients, up to the second harmonic as in (93), of 
the charge waveform giving the close lower bound to the maximum 
obtainable fundamental power, are 

a = l(* + p + t); b = 2d = £- (<r - r) ; 

(107) 
c = -2e = y^.(a + T - 2p). 

We will restrict ourselves to that class of voltage-charge relationships, 
V = V(Q), for which L(a,p,r) in (73) attains its maximum, subject 
to — m ^ (<r,P,r) ^ 1, when 

it = 1, r = -ra, -m < p < 1. (108) 

It would seem feasible that we might obtain a reasonable approxima- 
tion to the maximum power transfer from the fundamental to the second 
harmonic, by suitably shifting and expanding (or contracting) the 
above Fourier approximation, so that (94) is satisfied. Setting a = 1, 
T = — m in (107) and carrying out this procedure, we obtain the ap- 
proximating charge waveform 

Q(x) -i(l+ p -») +^?(1 +m)(2sinx +sin2z) 

(109) 
+ - (1 — to — 2p)(2 cos x — cos 2x) . 

If we define Q(x) = Q[x + (2x/3)], then (96) is satisfied and in the 
canonical representation of Q(x), (97), we have 

For the abrupt-junction diode operated in the region between forward 
conduction and reverse breakdown, p = ^, setting m = 0, v = 2 in 
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(100). Hence, from (110), s = -(1/V3) and y = 1/3, so that, from 
the previous section, the approximating charge waveform is actually 
the one which gives the maximum power transfer. 

5.4 The Numerical Computation of the Maximum Power Transfer, for 
Particular Diodes 

We have already obtained a two-parameter canonical representation 
of the charge waveform containing no higher than second harmonics 
and satisfying (94). The two parameters s and y lie in the bounded 
region given by (99), and Pi vanishes, independently of the voltage- 
charge relationship, on the boundary of this region. Also, since P x is 
antisymmetric in s, it is sufficient to consider only half the region and 
to maximize | Pi |. The maximization was carried out numerically for 
particular diodes, by means of the iterative process of fitting a quadric 
surface. As a starting point s (l) ,y (l) in the process, that point correspond- 
ing to the approximating charge waveform, derived in the previous 
section, was used. 

The results of the numerical computations for the voltage-charge 
relationship of (90), with v = 2 and v = f , and several values of m, 
are tabulated below. Tables II and III give the values of the maximum 
power transfer, max Pi , together with the corresponding values of the 
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max Pi 


Ri 


Ri 


/mux 


(* j + ^ 


(<* 2 + e^) 





0.2814 


1.462 


0.7310 


0.7698 


0.1482 


0.0370 


i 


0.7773 


1.966 


1.060 


1.160 


0.3289 


0.0865 
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1.549 


0.5947 


0.1573 
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1.451 


0.3484 


7 


3.366 


3.854 
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3.422 


3.616 


0.7414 


5 


4.371 


4.788 
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4.515 


6.869 


1.250 
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5.544 


6.020 
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20.95 


3.096 
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0.1623 
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0.7389 


0.7684 


0.1499 


0.0366 
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2.137 


1.182 
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0.3257 
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1.246 


2.428 


1.529 


1.560 


0.5635 


0.1652 
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2.271 


3.023 


1.882 
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0.3682 
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3.445 


3.479 


0.7703 


5 


4.691 


5.115 


1.907 


4.533 


6.710 


1.277 



704 



THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1962 



reactive powers in the fundamental and second harmonic, P : and P 2 , 
and the maximum current 7 ma x associated with the charge waveform 
Q(x), that is the maximum value of \Q'(x) |. It is worth noting that 
P 2 does not continue to increase with m. Also included are the squares 
of the amplitudes of the first and second harmonics in the charge wave- 
form, (b 2 + c) and (d 2 + e 2 ). These, together with the real and re- 
active powers, determine the normalized impedances. Tables IV and V 
give the values of — s and y which given max Pi , and also y ll) and Pi ' , 
the value of P, corresponding to y w and -s (1) = 1/V'S = 0.5774. It 
is interesting to note how close Pi (1) is to max Pi , except for the larger 
values of m. Table VI compares max P\ with the maximum obtainable 
fundamental power, max Pi , as obtained in Section III, for the case 
v = 2 and several values of m. It is also worth comparing the value of 
max Pi = 0.162 for the case v = $ , m = with the corresponding value 
of max Pi = 0.408. 



Table IV— (v = 2) 



m 


— X 


y(') 


y 


Pi<D 


max Pi 





0.5774 


0.3333 


0.3333 


0.2814 


0.2814 


i 


0.5839 


0.2963 


0.2942 


0.7770 


0.7773 


1 


0.5848 


0.2500 


0.2426 


1.283 


1.284 


2 


0.5716 


0.1852 


0.1782 


2.192 


2.198 


7 


0.5465 


0.1317 


0.1301 


3.307 


3.366 


5 


0.5246 


0.1019 


0.1046 


4.199 


4.371 


7 


0.5008 


0.0781 


0.0844 


5.197 


5.544 


9 


0.4810 


0.0633 


0.0717 


6.047 


6.586 







Table V- 


~{v = |) 






in 


— s 


y<» 


y 


;>i»> 


" max Pi 





0.5742 


0.3704 


0.3704 


0.1622 


0.1623 


j 


0.5871 


0.3566 


0.3562 


0.6775 


0.6782 


1 


0.5977 


0.2963 


0.2829 


1.241 


1.246 


2 


0.5875 


0.2112 


0.1989 


2.262 


2.271 


7 


0.5591 


0.1449 


0.1419 


3.518 


3.575 


5 


0.5331 


0.1097 


0.1130 


4.536 


4.691 







Table VI — (p = 


2) 






in 





i 


l 


2 


I 


5 


max P\ 
max P\ 


0.281 
0.687 


0.777 
1.83 


1.28 
3.02 


2.20 
5.50 


3.37 
9.33 


4.37 
13.15 
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5.5 On the Power Transfer From the Fundamental to the Third Harmonic 

Breitzer et al 3 were also concerned with the abrupt-junction diode 
operated in the region between forward conduction and reverse break- 
down and considered charge waveforms containing no higher than third 
harmonics. They treated in detail the power transfer from the funda- 
mental to the third harmonic, subject to P 2 = 0, and obtained a max- 
imum value of 

Pi = (0.0242)/ = 0.238 = -P 3 , (111) 

making allowance for the difference in notation. This value of Pi arose 
from two distinct charge waveforms. One was 

Q(x) = (0.5) + (0.310) sin a; + (0.168) sin 2x + (0.155) sin 3a;, (112) 

and the other was quite close to this. We saw previously how by taking 
the Fourier approximation, containing up to second harmonics, of a 
charge waveform which gives a good lower bound for max Pi subject 
only to restrictions on Q max and Q m m , and suitably shifting and ex- 
panding (or contracting) so that the restrictions on Q max and Q m j n are 
satisfied by the approximating charge waveform, we could obtain a 
good approximation to the maximum power transfer from the first to 
second harmonic, when no higher than second harmonics are allowed. 
In the case of the abrupt-junction diode operated in the region between 
forward conduction and reverse breakdown, which is the diode that 
we will consider in this section, it was found that the charge waveform 
so derived was precisely one that gives the maximum power transfer. 

Now, it is found that the best mean square approximation containing 
up to third harmonics, and subject to P-i = 0, to the charge waveform 
which gives the good lower bound to the maximum obtainable funda- 
mental power is 

Q(x) = fj + ^/(*)]> (H3) 

where 

f(a-) = [(0.4) sin .r + (0.25) sin 2x + (0.2) sin 3.rJ. (114) 
We shift and contract Q(x) by setting 



Q(x) = \ 



>+f 



M = max [/(a)], (115) 

so that Q I111IX = 1 and Q m \„ = 0. For this charge waveform, 
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P 1 = (0.0075)tt7M 3 = -P 3 ■ P 2 = 0. . (116) 

It is found that 

M = 0.680; Pi = 0.235, (117) 

and Q(x), as given by (114) and (115) is plotted in Fig. 7(a). The value 
of Pi in (117) is very close to the maximum value obtained by Breitzer 
et al, (111), and it is interesting to compare Fig. 7(a) with Fig. 7(b) 
which depicts Q(x) as given by (112). 
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Fig. 7 — Charge waveforms giving (a) approximately, and (b) exactly, the 
maximum power transfer from fundamental to third harmonic in abrupt-junction 
diode operated in the region between forward conduction and reverse breakdown. 
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VI. THE MAXIMIZATION OF THE POWER TRANSFER FROM THE FUNDA- 
MENTAL TO THE SECOND HARMONIC, FOR THE CURRENT-LIMITED 
DIODE 

6.1 The Canonical Representation of the Charge Waveform 

We are concerned with charge waveforms as in (93) and impose the 
restrictions 

Q ra .x=i, |Q'Ux = fc. (118) 

We observe that, for voltage-charge relationships of the form given by 

(90), 

maX [Pi | Qmax = P, | (/ | max = /] 

r /i < 119 ) 

= p ( " + " max \l\ | Q„ 111X = 1 , | Q Lax = -J, < p g 1 . 

In Appendix C we determine a canonical representation of Q{x), sub- 
ject to the conditions 

Qmax = 1; Q.nax ^ A'5 Qmin = ~fc, (120) 

by making use of the canonical representation obtained in Section 5.1, 
when the charge waveform has prescribed maximum and minimum 
values. Note that if Q(x) = Q(t - x), where Q(x) satisfies the con- 
ditions of (120), then 

frnax = 1; On- = ^ QLn ^ ~k. (121) 

From Appendix C, the five coefficients in (93) are given in terms of 
two parameters s and y. It is found that 

b = 7 ^-: 2d -r-^r; c=- r -*!L-4e, (122) 

(w - z) (w - z) (w — z) 

where 

M, = [i(l - S A ) - s 2 y}; z = \{y(l - s 2 ) - |(1 + s 2 ) 2 J, (123) 
and that 

a = [1 — max (b sin x + e cos x + d sin 2x + e cos 2a;)], (124) 

which in general has to be determined numerically. The waveform is 

translated so that 
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Q'(tt) = -k = QLin- (125) 

The parameter s arises from the equation 

Q'^tan- 1 s) = QLx ^ k, (126) 

and the parameter y is subject to the condition 

OSy ^ i(l - s 2 ), (127) 

which of course also implies that s 2 ^ 1. It is shown in Appendix C that 

P,s0 for */ = |(1 - s 2 ). (128) 

In order to maximize Pi subject to (118), it is sufficient, in view of 
the correspondence between Q(x) and Q(x) = Q(t — x) given by 
(120) and (121), to use the above canonical representation and to 
maximize | Pi |. 

(3.2 The Abrupt- J unction Diode 

We now consider the abrupt-junction diode operated in the region 
between forward conduction and reverse breakdown, for which the 
voltage- charge relationship is V(Q) = Q\ ^ Q fS 1. We first maxi- 
mize | Pi | subject to the conditions of (120), and suppose that A; is 
sufficiently small that Q m \ n ^ 0. Using the canonical representation ob- 
tained in the previous section, Pi may be expressed in terms of s and y. 
Omitting the details, it is found that | Pi | is maximized, subject to the 
restriction (127), for s 2 = $, y = 0. The charge waveform giving this 
maximum is 



whore 



Q(x) = 1 +k[S(x) - S lllux ), (129) 



^ ) = (4sin:r-sin2.,) i 
6 



It is readily verified that <S' limx = g = — <S IIlin , where 

-l^f-°- m (131) 

Thus Q min = ( I - 2gk), so that Q min ^ for 2gk ^ I. This Q(x) ac- 
tually gives a negative value of Pi , so that Q(x) = Q(ir — x) maxi- 
mizes P, , and it is found that 

max P, =%k 3 = 0.731 1A- 3 , for A- ^ J- = 0.681. (132) 
2l 2g 

Fig. 8 depicts S(t — x). 
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Fig. 8 — Shifted and normalized charge waveform for maximum power trans- 
fer from fundamental to second harmonic in current-limited abrupt-junction 
diode, with maximum current less than a critical value. 



The fundamental reactive power corresponding to max Pi is 
Rl =*£(!- gk)k >. 



(133) 



But, for the voltage-charge relationship V(Q) = Q 2 , the addition of a 
constant to the charge waveform does not affect Pi . Hence, if instead 
of requiring Q max = 1 we just require g Q(x) ^ 1, we have 



Rx = ^-ak 2 = 8.78ak 2 ; gk ^ a ^ (1 - gk). 



(134) 



6.3 The Optimum. Operating Frequency 

.So far, no discussion has been made of the angular frequency w of 
the actual periodicity of the charge waveform. We here consider this 
factor in the case of the abrupt-junction diode operated in the region 
between forward conduction and reverse breakdown. Now the physical 
limitation placed on the maximum current magnitude takes the form 



Q'(,r)\ t*-, 



(13r,) 



from (20). Also, the actual fundamental power pi is, from (9), propor- 
tional to coPi . We thus consider the maximization of wPi as oj varies, 
where the charge waveform Q(x), containing no higher than second 
harmonics, is subject to ^ Q(.v) ^ 1 and the condition in (135). We 
make use of results from Section V, as well as from the previous section. 
Thus, we define 

max [Pi | Q max = 1, Q min = -m] = P(m), (136) 
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and let the value of | Q' \ max for the charge waveform which gives P(m) 
be denoted by K(m). Then, remembering that the addition of a con- 
stant to the charge waveform does not affect Pi , since V(Q) = Q , we 
obtain from (95) and (103), 

P(m) =-i^=(l + m) 3 , (137) 

81'v 3 

and also, from (104), 

K{m) = (1 + m)K(O) = -^= (1 + m). (138) 

Similarly, we define 

max [P, | Q nmx = l,\Q' | ma x = k] = U(k), (139) 

and let the value of Qmin for charge waveforms which give U(k) be de- 
noted by — M(k). Then, from the previous section, 

H(Jb) = |£tf; M(k) = -(1 - 20b), (140) 

where g is given by (131). 

Now if Q(x) is subject to just the restriction ^ Q(x) ^ 1, then 
max P t = P(0), from (137). But, from (138), if («/«) ^ (3\/3)/4 
then the Q(x) which give this value of max Pi satisfy (135). Hence, 

w\ ,-, „,-»,-,,. /w\ _ . / CO 



^J max Pi = 0.2814 ( = 1 , g 1 = 1 ^ 1.299. (141) 

Note that if (w/k) > 1.299, then this gives an upper bound on 
(u/k) max Pi . Also, if («/«) > 1.299, then max P x ^ P(m) if tf(ro) = 
(*/«), and hence, from (137) and (138), 

^ max Pi ^ 0.617 f-Y , ^ ^ 1.299. (142) 



© 



From (140), setting fc = (k/ w )> we have 

(fj max Pi = 0.731 ^Y , (fj ^ 1.468, (143) 

and if ^ («A) < (1.468) = 2#, then this provides an upper bound 
on (u/k) max Pi . Also, if ^ (u/k) < 2g, then max P x ^ n[l/(2^)], 
from (140). Hence, 

-) max Pi ^ 0.231 Q , ^ Q ^ 1.468. (144) 
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Fig. 9 shows (w/k) max Pi as a function of (o>/k). For 



1.299 g I- J g 1.468, 



(145) 



the curve lies between the dashed lines. Thus, from the viewpoint of 
maximizing the actual fundamental real power, the optimum operating 
frequency, when the diode is not allowed to operate in the forward 
conduction region, lies in the range given by (145). Also, we can assert 
that 



0.3655 = — ^ max 



GH 



ol 



(146; 



0.4 Maximization of the Power Transfer, When the Region of Operation 
Includes Forward Conduction 

In a previous section we obtained a canonical representation of a 
charge waveform Q(x), containing no higher than second harmonics, 
for which Q MaK — 1, Qmax. ^ k and Qmia = — k. This canonical repre- 
sentation is given by (93), (122), (123) and (124), and involves two 
parameters s and y which lie in a bounded domain given by ^ y ^ 
5(1 — s )• It was shown that, independently of the voltage-charge 
relationship, Pi = on y = %(l — s 2 ). Moreover, it was seen that in 
order to maximize P, subject to Q Iliax = 1 and | Q' | max = k, it is suffi- 
cient to consider this canonical representation and to maximize | Pi |. 
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Fig. 9 — Maximum power transfer from fundamental to second harmonic in 
current-limited abrupt -junction diode operated in the region between forward 
conduction and reverse breakdown, vs. frequency. 
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The maximization was carried out analytically for the abrupt-junction 
diode when k is sufficiently small that the diode does not operate in the 
forward conduction region. We treat here, by means of numerical com- 
putation, the abrupt-junction diode when partial operation in the for- 
ward conduction region takes place, the normalized voltage-charge 
relationship being given by (61). 

Again, the maximization process was that of fitting a quadric surface, 
and this time it was also necessary to calculate a in (124) numerically. 
Further, it was desirable to first compute the value of Pi over a rough 
grid, and then to pick appropriate values s (1> and y w , as a starting point 
in the maximization process. Thus, for several values of k, max | Pi |, 
i.e., n(/v) in the notation of (139), was computed in the manner de- 
scribed above. For the values of s and y which gave max | Pi |, the cor- 
responding values of Pi and R 2 , the reactive powers in the fundamental 
and second harmonic, and of Q n ,in , i.e., — M(k) in the notation of the 
previous section, were calculated, together with (b + c ) and (d + e ), 
the squares of the amplitudes of the first and second harmonics in the 
charge waveform. The results of the numerical computations are tabu- 
lated in Table VII. We note that the values of Pi corresponding to the 
given values of s and y are negative. If Q(x) is the charge waveform 
corresponding to s and y, (93), (122), (123), and (124), then the posi- 
tive value of Pi , that is U(k), is obtained from the charge waveform 
Q(x) = Q(t — .r), or any translation thereof. 

Now, from (119) with v = 2, and from (139), 

max [Pi | Q max = p,\Q' Ux = I] = pH Q , < p ^ 1. (147) 
For Q max ^ we have Pi = 0, from (01). We may write 

p3n G) _ Q n W n(Q (148) 

n(A-) l~ s ll(l) U(k) " 

The quantity k~ 3 U(k) is depicted in Fig. 10(a), and it is seen to be a 
nonincreasing function of k. It follows, from (147) and (148), since 
U(k) is a strictly increasing function of k, that max Pi subject to Q max ^ 
1 and | Q' | ma x ^ k is attained with Q max = 1 and | Q' \ max = k. For 
A- < l/(2g() = 0.G81, it can also be attained with 2gk ^ Q max < 1 and 
| Q' | mux = k- We comment that for the voltage-charge relationship 
V(Q) = max (0,Q), max Pi subject to Q max ^ 1 and | Q' | max ^ k is 
not attained with Q max = 1, for sufficiently small k, since in this case 
Pi = if Q min ^ 0. 
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Let us now consider the frequency factor, as we did at the end of the 
previous section, so that (135) holds. Hence, setting k = (k/co). 



max U A J = max |_— _• 



(149) 



The curve in Fig. 10(b) depicts U(k)/k and it is seen to be an increas- 
ing function of k in the range shown, although it is to be expected that 
it tends to zero as k -*■ *> . It appears that max [U(k)/k] ~ 1, so that, 
from (146), a considerable improvement is obtained if the diode is per- 
mitted to operate in the forward conduction region. We must bear in 
mind, however, that we have idealized the voltage-charge relationship 
in the forward conduction region. 
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Fig. 10 — Maximum power transfer from fundamental to second harmonic 
divided by (a) the cube of the maximum current, and (b) the maximum current, 
for current -limited abrupt-junction diode with operation in forward conduction 
region permitted, vs. the maximum current. 
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Table VII 



k 


n(*) 


Ri 


Ri 


-»/(*) 


0.75 

1.0 

1.5 

2.25 

3.0 


0.3058 

0.6159 

1.265 

2.169 

2.979 


2.167 
2.440 
2.840 
3.483 
4.153 


0.3033 

0.5075 

0.8274 

1.058 

1.135 


0.0896 

0.3980 

1.027 

2.013 

3.024 


ft 


—s 


y 


(6 s + c») 


(</* + « 5 ) 


0.75 

1.0 

1.5 

2.25 

3.0 


0.5988 
0.6647 
0.7035 
0.7058 
0.6996 


0.0018 
0.0288 
0.0834 
0.1354 
0.1663 


0.2404 

0.3653 

0.7159 

1.613 

3.006 


0.0169 
0.0393 
0.1124 
0.2775 
0.5039 
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APPENDIX A 

From (28) and (29), for any X (which we take to be real), 



/', = 



(XQ(a:) - V[Q(x)]} sinxdx 



i\Q(,r) - V[Q(x)]\cosxdx 



Q(x) cos 



-a: 

= A [ [\Q(x) - Y[Q(x)]\ sin (x - 6) dx, 



s x dx J 
x) sin x dx) 



Q(x) s 



where 



(150) 



A sin 6 = I Q(x) smxdx; A cos 0= / Q(x) cos x dx. (151) 

Jo ■Ml 

Hence, 

A = / Q(x) cos (x - 6) dx. (152) 
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Now, max 1\ = max \P\\. Since 



/ f(x) sin (x — tp) dx 



^ 2 (max/ - min/), (153) 



(31), (150) and (152) lead to (85) and (86) in Section 4.2. We next 
derive the inequalities (88), where L is defined by (73) and (74). Now, 

L ^ max L(<r,l,—m) 

-m^irg 1 

= max [(1 + m)V(a) - (to + a)V(l) + (a - l)F(-m)] 

-'"^ 1 (154) 

h ■ x • MV(1) - V(-m)] 1r , J 

= -(1 + to) mm < / \ ^- J - F(<r) 

-m^Sl I (1 + TO) 

- [?nV(l) + 7(-m)]. 
Also, 
L ^ max L(l,p,—m) 

-m g p 5 1 

= (1 + „) max Uva)-y(- m) | _ 1 (155) 

-mgpgi ( (1 4- m) J 

4- [mV(l) + F(-to)]. 

Hence, from (86), (154) and (155), 

2L £ (1 4- to)£/. (156) 

Also, from (73) and (74), 

L = max ((t - p)[\<r - I'fcr)] 4- (<r - r)[Xp - V(p)\ 

-****** (157) 

+ (p - er)[Xr - r(r)]j, 

for any (real) A. In view of the remarks preceding (74) we may assume 
either that — to ^ a ^ p ^ t ^ 1, or that — to ^ t ^ p ^ <r ^ 1, 

without loss of generality. In the former case 

(r - p)[\a - V(c)] 4- (a - T)\\p - Y(p)\ 4- (p - <t)[\t - V(t)\ 
^ (t - p) max [\k — V{k)] 4- (o- - t) min [Xk - 7(<c)] 

+ (p - a) max [\k - V(k)\ (158) 

— m g * g 1 

= (t - <r){ max [Xk - I'U)] - min [X« - V(k)]}. 

— mfiiSl — m<ic<l 
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Hence 

(r - p)[Xa - V(a)] + (<r - r)[\p - V(p)] + (p - <r)[Xr - F(r)] 

g (1 + m){ max [X* - 7(k)1 - min k* ~ ^(«)11- (l 

Equation (159) may be derived, in a similar manner, when — m ^ 
r £ p £ * £ 1. Thus, from (157), 

Lg(l + m)| max [Xk - 7(«)] ~ min ft* - F(k)]}. (160) 

But this is true for all (real) X. Hence, from (8(5) 

L ^ (1 +m)tf. (161) 

If we do not restrict the voltage-charge relationship then the bounds 
given by (150) and (161) cannot be improved. This is demonstrated 
by considering the (somewhat artificial) relationship 

[ 1, Q= [(1 + m)a - m]\ 

V(Q)=<-1, Q = [1 - (1 + m)a]\ (162) 

0, otherwise; m > — 1, < a < J. 

It may be verified that in this case 

L = (1 + m) = (1 - «)(1 + m)f/. (163) 

We now find a class of voltage-charge relationships for which the 

bound in (161) is attained. If a ^ r, then, by the definition of U in 

(86), 

( a - T )U = ^ max [[7(a) - V(t)] P - (a - r)V(p)\ 

-** 1 (164) 

- min |[F(<7) - K(r)]p - (<r - r)K(p)}. 

— m g p g 1 

Let -m ^ r ^ <r ^ 1. Then, 

L ^ max [(p - t)F(» + (t - c)V( P ) + (a - P )V(t)\ 

—m <jp g I 

= max {[Via) - V(r)]p - (a - r)V(p)\ 

— mgpg 1 

+ [ctV(t) - rV{a)] 
£ (a - t)U + [<7F(r) - tV(v)] (165) 

+ min {[7((r) - V(t)) p - (a - r)7(p)j 

— m g p g I 

= (a - t)U 

+ min [(p - r)V(<r) + (r - <r)7(p) + (a - p)7(r)]. 

— mgpgl 
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Q(t/2) = [K3-2//)(l +m) -m)= [1 - 1(1 + 2/y)( 1 + m)J, (181) 

so that Q limx > 1 for s = 0, y < -$ and Q min < — m for 8 = 0, y > |. 

Hence Q nmx > 1 for y < 0, and Q lllin < —m for y > (I — s 2 ). The re- 
gion of interest, i.e., Q I1111X = 1 and Q min = —m, is given by 

g J 1(1 - s 2 ). (182) 

We next consider the fundamental power when the charge waveform 
Q(x) contains no higher than second harmonics. From (28), (29) and 
(170), 

Pi = * I V[Q(x)](b cos x - c sin x) dx. (183) 

Jo 

We determine conditions under which P x = 0, independently of the 
voltage-charge relationship V = V(Q). This is clearly the case if 
b = = c, or if Q(x), as given by (170), is a single-valued function 
of (b sin x + c cos .r), for then the integrand in (183) is the derivative 
of a periodic function. Noting that 

2(6 sin x + c cos.r) 2 = (ft 2 + c 2 ) + 2bc sin 2x + (c - b' 2 ) cos 2x, ( 184) 
it follows from (170) that the latter condition holds if 

d = 2X/>c; e = X(c 2 - b 2 ), (185) 

for some X. Combining this condition with b = = c, 

2bcc + d(b 2 - c 2 ) = 0=> Pj = 0. (186) 

We now consider the canonical representation of Q(.r), with Q nmx = 1 
and Qmin = — m, wherein the coefficients in (170) are given by (97). 
Then condition (186) becomes, upon reduction, 

y = 0, or y = (1 - s 2 ), or s = =* P t ss 0. (187) 

APPENDIX C 

We here determine the canonical form of Q(x), as given by (170), 
such that 

Qmax = 1; Qmax ^ /."I Qmin = Q'(t) = "A". (188) 

_ Now, when the charge waveform Q(x) is subject to Q max = 1 and 
Qmin = —m, the five coefficients corresponding to those in (170) are 
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given in terms of two parameters s and y, from (97), by 

a = [(c - e) — m]; b = 2d = (1 + m)sy; 
c = (l + m)w; e = (1 -f m)z, 

where 

w = [i(i - s 4 ) - 8 2 y]; z = i[/y(l - s 2 ) - }(] + s 2 ) 2 ]. (190) 

The charge waveform is translated so that Q(ir) = —m, which may be 
done without loss of generality. The parameter s arises from the condi- 
tion Q(2 tan -1 s) = 1, and the parameter y is subject to the condition 
0^*/^(l—s"), which of course also implies s 2 ^ 1. If, in addition, 
a = 0, then 

(1 + m)(w - z) = m, (191) 

and hence, from (190), 

2y(\ + 3fi 2 ) = [(1 + s 2 )(5 - 3s 2 ) - 8m/ (1 + m)]. (192) 

Now ^ y ^ ( 1 — s 2 ) , but if we require ?n ^ 1 then 

g y £ i(] - A ( m ^ 0- (193) 

Turning to a charge waveform Q{x), as given by (170), which satis- 
fies the conditions of (188), we may write 

Q'(x) = -Q(x), mfe 1, (194) 

m 

where (191) and (193) hold. Hence, 

,, v _ k[s y( sin x + \ sin 2x) + w cos x + z cos 2.r ] fig.-) 
* ' " (w — z) 

Integrating, and remembering that Q max = 1 , 

Q(x) = (I +k[S(x) - S, nax \\, (19(5) 

where 



\w sin x + - sin 2.r — sy(cos x + J cos 2z) (107) 

S(x) = - t : . 

(w - z) 

In general, «S max = max[jS(.r)] is determined numerically. 
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We now turn to the fundamental power, P x , when Q(.v) has the above 
canonical representation. From (170), (186), (190), (196) and (197), 
we find that P x = 0, independently of the voltage-charge relationship, 
if 

[(1 - s 2 ) - 2//]{(l - s';)(l + 3s 8 ) - s\{\ - i) - 2yf\ = 0. (198) 

In view of (193), the second factor vanishes only if s"' = 1, y =*0 
Hence we conclude that 

P, =~ for // = !(!- s 2 ). (199) 



REFERENCES 

1. Early, J. M., Maximum Rapidly-Switchable Power Density in Junction 
Triodes, IRE Trans, on Electron Devices, ED-6, July, 1959, pp. 322-325. 

2. Early, J. M., private communication. 

3. Breitzer, D. I., Gardner, R., Greene, J. C, Lombardo, P. P., Salzberg, B., 
and Seigel, K., Third Quarterly Progress Report, Application of Semicon- 
ductor Diodes to Low-Noise Amplifiers, Harmonic Generators, and Fast-Acting, 
TR Switches, Airborne Instruments Laboratory. Report No. 4589-1-3 (March, 
1959). 



